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2'he Thermal Stresses in Spherical Shells Concentrically Heated.. 

By CriAiiLES H. Lees, D.Sc, F.R.S. 

(Received June 8, 1921.) 

lNTPvOr)UCT10xN[. 

1. In recent years there have been a number of failures of materials when 
unequally heated which appeared to be due to the stresses set up by the 
differences of temperature in the materials. These stresses, although 
unimportant for small differences of temperature and at low temperatures, 
may be serious for large differences at high temperatures, owing to their 
increase with, the temperature differences and to the great reduction of the 
elasticities and ultimate strengths of materials at high temperatures. It is 
therefore desirable, when possible, to calculate the thermal stresses to which 
structures are liable to be subjected by the differences of temperature 
likely to exist in them, especially if these differences are large and the 
temperatures themselves high. 

The circumstances producing high thermal stresses are present to an 
exceptional degree in metallurgical, pottery and glass furnaces, and it is in 
these cases that disintegration of the materials used in the construction of 
the furnaces is most marked. The furnaces themselves vary so much in 
shape and material that it is advisable at present to discuss only some 
simple form and material, the investigation of which presents the minimum 
number of mathematical difficulties. 

The roughly hemispherical form of many furnaces suggests the considera- 
tion of the comparatively simple case of a concentric spherical shell 
concentrically heated, composed of materials for which the physical 
constants over the ranges of temperature considered may be taken as 
approximately constant. 

The case of a solid sphere concentrically heated was treated by Hopkinson*' 
in 1879 and by Almansif in 1897, on the assumption that, throughout the 
range of temperature considered, the expansion of the heated material was 
proportional to the rise of temperature. 

It is, however, possible to cover those cases in which the differences of 
temperature are too great to allow this assumption to be made, by dealing 
in the equations connecting strain and stress directly with the expansion of 
the material from some standard temperature to the actual temperature, 

■^ Hopkinson, J., 'Mess, of Maths.,' vol. 8, p. 168 (1879). 
t Ahuansi, ' Atti Accad. Torino/ vol. 32, p. 701 (1897). 
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instead of introducing the temperature and the coefficient of expansion. 
In the following deduction this more general method is adopted. 



Expansion of the Shell. 

2. Let r (fig. 1) be the distance of a point from the centre of curvature of a 
concentric spherical shell at a uniform temperature and subjected to no 
stress, and let the temperature at all points of a thin shell of radius r be 

raised by 0, where ^ is a function of r, 
and the pressure at the inside surface 
of the shell of radius ri become Ei, and 
that at the outside surface r2 = Ks. 

If the law lof linear expansion of the 
material of the shell is expressed by 
the formula Iq = lo{l-\'aO-j-'b6^-i-c0^, 
etc.), where a, h, c are constants, we 
may write ^ = ^o (1 + %), where 
ae = ad + hd^ -h c6'^ + etc., and % is then 
the actual expansion in length of a rod 
initially 1 cm, long when its tem- 
perature is raised by 6. In the same 
way, the change of volume may be 
represented by the formula Vs = % (1 + a^), where a^ is the dilatation of 1 c.c, 
when its temperature is raised by 6^ C. 

If a thin concentric shell of internal radius r and external r + dr he 
removed from the main shell while at its initial temperature, have its tem- 
perature raised so that the dilatation of 1 c.c. is a^, its volume, 4:7rrMr, will 
be increased by 4:7Tr^ocf^dr, 

If each elementary shell of the material be isolated in this way and suffer 
dilatation at its own special rate, the whole increase in volume of the 

elementary shells will be ^kwr'^ocs clr and the mean dilatation of the 




Fig. l.—The Shell. 



''■'>. 



4:7rT'^aQ dr 



material of the shell per c.c. *^''i 

In general, the expanded shells would no longer fit together, and it is in 
forcing a fit that the stresses known as thermal stresses arise. 

We note, in the first instance, that if the dilatation per cubic centimetre 
were the same for each shell, they would still fit together, and the thermal 
stresses would be zero. We therefore decompose the dilatation per cubic 
centimetre, a^ of each elementary shell in the general case into two parts : 
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first, a part, a, equal to the mean dilatation for the whole shell ; second, a 
part, a/, equal to the excess of the actual dilatation over that mean. The 
first part of the dilatation introduces no stresses ; the thermal stresses are all 
due to the second part.* 

The Stresses in the Shell. 

3. To find the nature and magnitude of the stresses and strains in the shell. 
Consider a small portion of the shell bounded by concentric spherical surfaces 
of radii r and o''-\-dT respectively, and by the surface of a cone of small vertical 
angle with its apex at the centre of curvature of the shell (fig. 2). Let E be 




R+(^R 



T 



T 

Fig. 2. — The pressures on an elementary disc. 

the radial pressure on the element of the surface of the disc having radius r, 
and E + c^E that on the surface T-\-ch\ Let T be the tangential pressure on 
the conical surface of the disc. These are by symmetry the only stresses. 
Let r' be the radial displacement outwards from the centre of curvature of 
any point initially at distance r from that centre. Again by symmetry this 
is the only displacement. 

The stress-strain equations for the elementary disc are : — 

(l_^)T-crE = ~6(~-aA (3.1) 

E-~2o-T= -e^^-aA (3.2) 

where e is Young's Modulus, a Poisson's ratio, supposed to be constants, and 
ae the excess of the linear expansion of the material of the disc over the 
mean expansion for the whole of the shell. 

■^ The absence of this limitation renders the equations of Rayleigh, ' Phil. Mag.,' 
voL 1, p. 170, No. 3 et seq. (1901) ; and of Love, 'Theory of Elasticity,' 2nd ed., p. 106, 
No. 34, too general. 
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The equation for equilibrium of the forces acting on the disc is 

2Tr =: -f- (Er^). 
dr 



(3.3) 



Substituting the value of T from (3.3) in the former equations they 
hecome 

(3.4) 



( 1 — cr) -— . -r- (Er^) — <tE 
^ -^ 2r dr ^ 



r , 

6 I % 

r 



and 



r dr 



clr^ 



\ dr 



-^— -a 



e 



(3.5) 



Multiplyiu^ (3.4) through byr, differentiating with respect to r and adding 
the result to (3.5) we eliminate r' and obtain on reduction 






(%€(/a 



dr 



. r 



dr 1 — cr . d/t 



ithe solution of which is 



B' 26 1 r 

E = A4--~Tr+ :i — - — q rhi0 . 6lr, 
r"" 1 — cr 7-"^ 



(3.6) 



(3.7) 



where A and B' are constants whose values depend on the pressures exerted 
on the inner and outer surfaces of the shell 

This expression for E suggests the substitution of v the volume of the 
shell of radius r instead of r in the equation, which becomes 



A B 2 € 

E.=: A+-4-K 



V 9 1 



o- V 



1 p 

a0 ,_ dv, 



•or 



'Rv = A'?; + B + 



2 



6 



Cv 



9 1~-cr 



a0' dv. 



(3.8) 



If El is the pressure on the inner surface n of volume Vi, and E2 that on 
the outer rg of volume % we have 



2 € 

El?/'! = Ai^i + B + K" 



'», 



and 



9 1-a 

2 6 

Eg^ = A% + B + -^ 



fif^ 



' . <2v 



9 1 



cr 






On eliminating A and B from these equations and reducing we have 



E =: E2+ (Ei-^Es) 



._, . %/'?;— 1 2 6 



t^g/'?^]— -1 9 ' 1 — (J ' -y 



■% 



(5e^^ cl'y, 



V 



"Jjf 



•Vr, 



writing E' for the expansion a/ dv ^ — \ a/dv, due to the rise of 



temperature of that part of the shell within the spherical surface of volume v 
we have 

E = E2 + (Ex — Eg)— ^1 r + ;^^r -.2—, (3.9) 

V2IV1 — I 9 l"-cr ti 



Sjyherical Shells Concentrically Heated, 383 

or if «t = Ijv 

E = ^^^~-^+ "li-Es) ™-^^+ J T-^ . 2^^E^ (3.90 

Ui—%i2 ^ Ui—'ifj2 9 l — cr 

To deterinine T we have from equation (3.3) T = E+- —--, and on 



Zi (h') 



introducing" the volume i? instead of the radius r, 



T==E + 2t?^. (3.10) 

2 d'V 

On substitutiug the values of E an^^ dR/dv in this equation and reducing 



we get 



T = E.-(E,-:E.)!^E/^ + J-l-(3./-n (3.11) 

V2/V1—I .9 1 — cr\ v; ^ ^ 

or in terms of zi 

E2'?^l — Ih'^2 In — E2 '2^ 1 6 tij ' .,Tr'/\ /ot-i/\ 

== ^ 4- ^^™_ (^^ijf^ —^6jii ). (3.11 ) 

tl^l — 2^2 ^ %i— -1^2 9 1— GT / X / 

The "mean pressure" or "equivalent hydrostatic pressure" to which the 
material of the shell at any point is subjected is (E4- 2T)/3 and the preceding 
equations give us 

I^+2T E2%— •Ei2^2 , 1 6 o / /OiON 

_— — — J^ — ^ . 2a^. (312) 

The two shear stresses in diametral planes at right angles to each other 
are (E— T)/3 and our equations give 

Writing t^' for the increment of volume 1; of the spherical surface of 
radius r owing to the displacement t\ and %' for the increment of '^ = Ijv, 

we have since ^ =™^ = — !L on substituting the values of E and T 
in equation (3.1) and rearranging 

r 3 -y 3 w 6 ' 'Zii— 2/^2 



It may be noted that 



l + <r El — E2 W ,11 + 0- E /oi.x 

-^ , ,^ _ ^ (0.14:) 

e 2i %— ^e2 9 1 — 0- -i? 



X — 2cr 1 , l + cr 1 

— and 



where 7i; is the bulk modulus and /a the rigidity of the material. 

The connection between r'/^' and E is obtained readily by eliminating T 

and as' from the original equations which gives ~ (^ = iif — . (3.1 5) 

* dr\r 2e dr ^ ■ / 
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The stresses and displaeemeBts due fco the internal and external pressures 
El and Eg agree with Lame's results.* 

On putting Ti = 0, Ei = E2 = and a^' = 3ad^ in these equations, the 
thermal stresses reduce to Hopkinson's values for a solid sphere in which the 
range of temperature is small enough to allow the expansion with rise of 
temperature to be taken as uniform. 

4 From these expressions it is seen that the cubical compression due to the 
pressures Ei and Eg on the surfaces of the shell is tlie same throughout the 
material,f while that due to the temperature distribution is at any point 
proportional to the excess of the thermal dilatation at that point over the 
mean thermal dilatation of the whole material of the shell. 

The radial pressure E, the displacement -r/i% the tangential pressure T, and 
the shear stresses in diametral planes at right angles to each other, due to the 
pressures on the surfaces, are all linear functions of it the reciprocal of the 
volume D of the sphere through the point. The radial pressure E, and the 
displacement //r, due to the temperature distribution are proportional to the 
quotient of the excess cubical expansion of the elementary shells between the 
inner surface th and the surface v, by the volume v of the latter surface. The 
other thermal stresses are found by the usual combinations of radial pressure 
and the mean pressure. 

We see from, the foregoing investigation that the general problem of the 
stresses in a spherical shell, produced by change of temperature such that 
isothermal surfaces are concentric spheres and by pressures applied to its 
surfaces, resolves itself into the determination by the methods indicated above 
of the stresses and displacements due to the expansion «/ of each elementary 
shell above or below the mean a for the whole shell, and the superposition on 
these of the displacements due to the mean expansion U and of the stresses 
and displacements due to the pressures Ei and E3 applied to the inner and 
outer surfaces of the expanded shell 

Graphical Eipbesentatiok of Steissis. 

5. It is not easy to realise the meanings of the foregoing expressions for 
the stresses due to mechanical and to thermal causes without graphic repre- 
sentation and the expressions themselves suggest the following method 
of representing the stresses at a point, r, in terms of the volume, v, of the 
sphere through the point, or in the first case in terms of the reciprocal, % of 
that volume. 

* See e.g.^ Love, ^ Theory of Elasticity,* 2nd ed., pp. 139, 140 (1906). 
t The constancy of the cubic compression due to the pressures on the two surfaces 
is noted by Lamh in his * Statics,' p, 332 (1912). 
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A. Stresses due to Pressures on Inside and Outside Surfaces, 

Set off along the axis of abscissae (fig. 3) the reciprocals OUi = ui and 
OU == U2 of the volumes Vi and V2 



So 



So 



}r^^^'^ 


^^^-"-"^^-^^ i 




"^^---^1 










of inside and outside surfaces. At 

Us erect an ordinate, ^<^% equal 

on some convenient scale to the 

pressure, E2, on the outside surface. 

Join the top, E2, of this ordinate to 

the point Ui and produce to cut the 

pressure axis in So. Through So 

draw a line parallel to the axis of 

abscissae, and let it cut the ordinates 

through Ui and U2 in Si and S2. 

Set off above Si a length, SiRi, equal to Ei, the pressure on the inside 

surface, and join the end, Ei, of this ordinate to S2, and produce to cut the 

axis of pressures in So'. 

Then OSo is E2 — ^^ , So So' is Ei — 5M_ : hence 



Fig. 3. — The mean pressure OSq' determined 
by U2R2> U1E2S0, SqSi, SjBg, and RiSgSo'. 



i I 



OSo 



E2t^l — E2t^l 



(5.1) 



that is the ''mean pressure" or "mean hydrostatic pressure *' throughout the 
shell. 



To get the terms involving (Ei — E2) 



u 



necessary to express the 



U1--U2 
various stresses, erect at Ui, fig. 4, the ordinate UiA = (Ei"~E2), join the 



SB 




U Vj 



Fig. 4. — UiA = Rj ~ Eg, OB is parallel to U2A1, US gives the stress at U. 



extremity A of it to the point U2 and draw through. .0. a straight line OB, 
parallel to this join, cutting UiA in B. The ordinate US at any point U is 



then (El — E2) . 



to 



. Hence the sheer stress in each of two mutually 



III — "^h 
perpendicular diametral planes through the point U 



1 Tiq 



(5.2) 



VOL. C. — A. 



2 c 



386 



Prof. O. H. Lees. The Thermal Stresses in 



For the other stresses we have 



Xt = + ( hi — ±t2 ) ~ 

U1—U2 '^^'l— '?^2 



m ^2'^! — 3ai^2 Rl — R2 

^/^ /2/_ 9 



v. 



lll—U2 

and for the displacements 



III — ^^2 



OSo' + US, 



OSo'--iUS, 



It 

3 V 



1 — - 2 cr Es'^^i — Ei?i2 , 1 + o- El — Eg %i- 



~\ 



6 



Ul — U2 

1 + cr 



2 '2(1 — '^/.2 



r' 



Obo 4' 



JUS 



(5.3) 



(5.4) 



(5.5) 



^ 



fc>' 



B; Stresses due to the Deviation of the Dilatation due to Tenvperatute from the 

mean Dilatation. 

6. Set off (fig. 5) as abscissa OVi, OV2, the volumes, t^i and 1)2 of the inner 
and outer surfaces of the shell Take as ordmates the values of E', the 
increase of volume of the series of infinitely thin shells between vi and v due 
to the excess of the temperature of each shell over the mean temperature of 
the whole of the shells between the inner and outer surfaces. 

We thus ^et the curve of volume increase, Vi, PY2, extending from Vi, 

where it cuts the axis of volumes, to Y2, 
where it again cuts the axis. As we 
proceed outwards it slopes upwards, 
where the layers are hotter, or the 
dilatations greater, than the mean, and 
downwards, where they are colder, or 
where the dilatation is below the mean. 
The slope upwards at any point is pro- 
portional to the excess dilatation per 
cubic centimetre on the scale taken. 
Join P, a point on the curve, " to the 
origin, 0. Through draw OQs parallel to the tangent to the E^ curve at P, 
and at V^, corresponding to a volume Vg, a convenient multiple of 10, erect 
a perpendicular to the volume axis, cutting OQs in Q^ and OP in P^. 
Then the radial pressure E at the point P is 




Fig. 5. — E' the excess expansion curve 
gives at P with OQ^ parallel to the 
tangent at P, the lines V^Ps, VsQs 
determining the stresses. 



9 1 



2> 



V.P 



s*- s 



V, 



(6.1) 



the " mean pressure " at P is 



V.Q. 



9 1 



cr 



V, 



(6.2) 
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the tangential pressure T is 



the two shear stresses are 



1 e 3VA~V,P, 



9 1 



(T 






9 1— o- Vg 
and the radial displacement is given by 

r' _lv' _1 1 + a YsVj 



(6.3) 



(6.4) 



S V 9 1 



(6.5) 



To the latter must be added in case the mean temperature has been raised 6 
the radial displacement given by — = ^ = - ^ due to it. 

Construction of the Excess Expansion Curve E'. 

7. The problem of determining the thermal stresses at a point in the 
material of the shell thus reduces to the construction of the curve E', and 
the drawing of the radius vector and tangent to the curve at the point. 

To draw the curve most conveniently we have E' = I a/ c^-y, where a/ is the 

excess of the cubical dilatation of 1 c.c. at 6, over the mean dilatation a 
for the whole shell. 

Hence if ae is the actual dilatation at 6, 



E 



rv Cv 

' = 1 (a -^oi^dv = aQdv — dLiy — Vi), 



or if 



E 



r 



dv 



and 



E, 






dv = 5t('y2— •'Vi), 



== Jtli — Jii2 . 

V2 + Vx 



(7.1) 



We, therefore, proceed as follows :- — Draw first the curve of temperature 
distribution, 6 (fig. 6), throughout the shell as determined by the temperatures 
^t the two surfaces and the laws of conductivity. 

D 




L 







\ 




Fia. 6. Fig. 7. 

Fia. 6. — B the curve of temperature, D the curve of dilatation. 
Fig. 7. — The excess expansion curve E' obtained from the expansion curve E. 

2 2 
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For each temperature estimate from tables of the expansion of the material 
of the shell the dilatation, uq, of 1 c.c. at each part of the shell, and draw the 
curve of dilatation per cubic centimetre, D, to correspond to the temperature 
curve. Find the total expansion, E, that is the area included between the 
axis of volume, the D curve, the ordinate through vi and that through v 
for different values of v up to %, and plot the values in terms of v^ curve E 
(fig.7). 

Join the end of the curve at V2 with its starting point at Vi. The excess of 
the ordinate of the curve at any point over the ordinate of the join is the 
ordinate of the excess cubical expansion curve, E', required. 

Deductions feom the Excess Expansion Curve. 

8. The curves of excess cubical expansion, E' in terms of v allow the 
conditions for zero and maximum values of the various stresses to be seen by 
inspection in many cases. If the dilatation due to temperature decreases as 
we proceed outwards through, the shell, the curve E' will have one maximum 
between its zeros at ^i and V2, while if the temperature and dilatation increase 
outwards E' will have one minimum between its zeros. If the dilatation is a 
maximum at some layer within the shell the curve, E', on starting from its 
zero at vi, is first negative, then passes through zero and becomes positive and 
ends at zero at V2. If the dilatation is a minimum at a layer within the shell, 
the curve starts by being positive near its Vi zero, passes through zero and 
becomes negative before reaching its V2 zero. 

In all cases the maximum and minimum ordinates of the E' curve occur at 
the layers for which the dilatation due to temperature has the mean value 
for the whole shell. The ordinate is a maximum if the dilatation decreases 
outwards and a minimum if it increases. 

The displacement is zero for points at which the E' curve crosses the axis 
of volumes. Two of these points are at the inner and outer surfaces 
respectively. The maximum outwards displacement occurs at the layer for 
which the tangent to the E^ curve where it is positive passes through the 
origin, and the maximum inwards displacement at the layer for which the 
tangent to the cvirve where it is negative passes through the origin. Both 
layers are nearer to the centre of the sphere than the layers of maximum or 
minimum E'. 

The radial pressure is also zero at points at which the E' curve crosses the 
axis. Two of the points are at the inner and outer surfaces respectively- 
The maximum value of the radial pressure occurs in the same layer as the 
maximum outward displacement, and the maximum radial tension in the 
layer of maximum inward displacement. These layers are a little nearer the 
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centre of the sphere than those of maximum and minimum ordinates of the 
E' curve respectively. 

The mean pressure about a point vanishes at the layers for which the 
dilatation is the mean for the whole shell, that is at the maximum E'. 
The maximum valu.es of the mean pressure occur in the layers for which the 
temperature dilatation differs most from the mean, that is, where the slope of 
the E' curve is greatest. In most practical cases this will be at or near the 
inner and outer surfaces and at the layer within at maximum or minimum 
temperature if such exists. The mean pressure is positive in the hotter and 
negative {i,e,, a tension) in the cooler portions of the shell. 

The two shears in diametral planes at right angles to each other are zero for 
the layers for which the tangent to the E' curve passes through the origin, 
that is for the layers for which the displacements and the radial pressures 
are numerically greatest. They are greatest when the temperature dilatation 
is above the mean and is increasing outwards, or when it is below the mean 
and is decreasing outwards. 

The tangential pressure is zero in a layer between that of mean dilatation 
and that of maximum radial pressure, for which the slope of the E' curve is 
one-third that of the join of the point on the curve with the origin. In 
general the tangential pressure is positive and large at or near the layer at 
highest temperature, and negative {i.e.^ a tension) and large at or near that at 
lowest temperature. 

From the foregoing it is seen that the determination of the stresses produced 
in a spherical shell by differences of temperature concentric with the shell is 
reduced to that of calculating the expansions of the elementary layers of the 
shell due to the deviation of their temperature from the mean temperature of 
the whole shell. 

The Temperatuke Distkibutiok. 

9. The determination of the distribution of temperature throughout the 
shell, given the temperatures of the two surfaces, is a problem in heat or 
temperature conductivity, which if the conductivity of the material were 
independent of the temperature, could be solved by the usual methods. If 
the thermal conductivity, h, of the material cannot be taken as independent 
of the temperature the condition for continuity of flow of heat in the material 
must be given the general form 

where c is the specific heat and p the density of the material. 

The difficulty of determining ^ as a function of the time and radius by 
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means of this equation, when c and h both vary considerably with temperature, 
is sufficient justification for making some simplifying assumptions in order to 
obtain some idea of the magnitudes of the thermal stresses in materials under 
common practical conditions. 

The Steady State. 
10. In the first instance let us assume that the temperature at each point 

of the material has become steady. Our equation reduces to — - ( /^r^ — ] = 0, 
which on integration gives 

dr 4:iTr^ 

where H is the heat which flows outward through the shell per second. 

The heat conductivity k in general decreases with increase of temperature,* 
but there appear to be some furnace materials for which it is approximately 
independent of temperature.f 

Assuming that the shell is of such a material, we have for the temperature 

If the volume v^—vi of the material of the shell is small compared to its 
internal volume vi this may be given the approximate form 

0-0 2 "02— V /in9\ 

01 — 02 ^2 — ^1 

where as before, 0i and 02 are the temperatures at the surfaces and the 
mean temperature ■ 



^1 + 02 



2 

Taking the expansion of the material of the shell to be proportional to the 
excess of temperature over the mean, this gives for the expansion, E', of the 
elementary shells between the inner surface, vi, and that of volume v, 

The maximum value of E' is -' (^i — ^2) ' ^'^"T"^^ and occurs at 'y = -— -— i. 
The maximum R is approximately 

i-^.^(^i-^.)^^^^, (10.4) 

9 1 — <r 2 V2-i-Vi 

^ Lees, ' Phil. Trans. Eoy. Soc./ A, vol. 191, p. 399 (1898) ; Eucken, 'Ann. der Phys.; 
vol. 34, pp. 217, 219 (1911). 

t Clement and Egy, ' Univ. of Illinois Bulletin,' vol. 6, No. 42 (1909). 
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and occurs slightly within the middle layer of the shell. It is a pressure if 
^1 > ^2 and a tension if 6i < 62- 
The maximum T is approximately 



1 - ^^ 



O <5t 



and is found at the surfaces. At the inside it is a pressure, at the outside 
surface a tension if di > 62- 

The maximum value of the mean pressure about a point is 

\^.u{e^-e,), (10.6) 

and occurs at the surfaces. 
The maximum shears are 



at the surfaces. 

The maximum displacement is 

11 + '- «(^,_^,)H?^, (10.8) 

and occurs approximately at the middle layer of the shell. 

NuMEKicAL Values : Practical Steady Case. 

11. When we attempt to assign numerical values to these quantities in 
any practical case, we are met by the difficulty that we have little or no 
knowledge of the expansion of the materials of furnaces at high tem- 
perature, and no knowledge whatever of their elastic constants and ultimate 
strengths at high temperatures. 

If the wall of the furnace is 20 cm. thick, and the temperature inside is 
1100° C, and outside 100° C. 6^1-^2 is 1000° C, and, if we take the wall to 
be of material of eibout the physical properties of jSrebrick at ordinary 
temperatures, we should have the coefficient of cubical expansion 12 x 10"^, 
Young's modulus about 3 x 10^^ Poisson's ratio about 1/3, a crushing 
strength of about 5 x 10^ and a tensile strength of about 2 x 10^ dynes 
per square centimetre.* At higher temperatures the expansion is greater, 
Young's modulus less, the crushing strength less, and the tensile strength 
probably considerably less than at ordinary temperatures. 

Calculating the stresses on the basis of the above figures we have for the 
maximum : — 

Eadial pressure E = 3 x 10^ . ^^"~^^ near the middle of the wall. 

V2 + Vi 

* From Mark's * Mechanical Engineer's Handbook.' 
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Tangential pressure T == 3 x 10^ positive at the inside surface and 
negative, that is a tension, at the outside surface. 

Mean pressure at a point = 2 x 10^, positive at the inside, negative at the 
outside surface. 

Shear stress = 1*0 x 10^ at the surfaces. 

Displacement = 0*7 x 10"^ ^^^;Z^. near the middle of the wall. 

These maxima exceed the limits of tensile strength of the material, and 
cracks would be formed on the outer surface, even if the wall were solid. If 
the bricks were set in cement, the cracks would form at the joints. 

Unsteady States. 

12. If the temperature throughout the shell is not steady, the W curves 
are much more complicated in shape, and some of the stresses may be 
considerably greater than they are in the steady state. 

Suppose, for example, that the shell is all at 100° C. and the inner surface 
is suddenly raised to 1100° 0. After the first few seconds, the E' curve, 
corresponding to the temperature distribution, starts from the vi zero at an 
angle to the axis of volume whose tangent is double that for the steady state. 
The maximum ordinate of the curve is, however, much smaller than in the 
steady state, and the inclination of the curve to the axis of volume at the V2 
zero is also much smaller. Hence the tangential pressure, T, the shear 
stresses, and the mean pressure at the inner surface, are approximately 
double those in the steady state at the inner surface, and the tangential 
tension at the outer surface smaller than that in the steady state. The 
radial pressure, E, and the radial displacement do not exceed the values 
attained in the steady state. 

As a second case, suppose the shell, after attaining the steady state with 
inside temperature 1100° C. and outside 100° C, has its inside surface 
suddenly cooled to 100° C. The E' curve at a second or two after the 
change proceeds from the inner Vi zero, with a downward slope equal to the 
upward slope in the steady state. The tangential pressure at the inner 
surface in the steady state becomes an equal tangential tension, the mean 
pressure and the shear stresses are reversed, while the maximum radial 
pressure and the displacement retain approximately the values they had in 
the steady state. At the outer surface the stresses change from their steady 
values slowly. 

Thus a sudden rise of the temperature of the inside of the shell from its 
initial value, supposed to be that of the outside surface, introduces at the 
inner surface double the tangential pressure, the mean pressure and the 



Spherical Shells Concentrically Heated. 393 

«hear stresses which exist at that surface in the steady state, hence spalling 
may result. 

A sudden fall of the inside temperature from that of the steady state 
to that of the outside surface converts the tangential pressure of the steady 
state at that surface into an equal tension, the mean pressure into an equal 
tension, and reverses the sign of the shear stress. 

Although the tangential pressure at the ner surface in the steady state 
may not have been sufficient to cause crushing, the tension thus produced 
may cause radial cracks, owing to the tensile strength of the material being 
less than its crushing strength. 

Furnaces. 

13. Through the kindness of Sir Eobert Hadfield and of Dr. J, W. Mellor 
I have been able to inspect furnaces which have been in use for some 
time and to examine fractured materials, and the appearances of the cracks 
in the bricks or other materials of the furnaces show that the above theorv 
reproduces the principal facts. The general tendency of rapid internal 
heating to produce spalling of the inner surface, and of rapid cooling to 
produce radial cracks, is quite well known to those in charge of the furnaces, 
and proper steps are taken to reduce the rate of change of temperature 
when the furnace has to be heated or cooled. A furnace which is to be 
taken down or repaired is cooled rapidly by blowing cold air through it in 
order to produce disintegration of the walls and facilitate the taking-down 
process. 

Conclusion. 

14. It is thus evident that the stresses in the wall of a furnace of spherical 
or approximately spherical form produced by a distribution of temperature 
throughout the wall, which is either known from observation or calculable 
hy the laws of conduction of heat, can be determined readily by graphical 
^or other methods if the law of expansion of the material throughout the 
a^ange of temperature is known, and Young's modulus and Poisson's ratio 
or any two elastic moduli are also known, and can be taken as approximately 
•constant. Whether the stresses thus determined will produce crushing 
or rupture of the wall depends on the values of the crushing load or the 
tensile strength of the material at the temperature which exists in the 
region of the stress. But the forms of the expressions for the stresses show 
that these stresses are proportional to the product of the dilatation of the 
material with rise of temperature into Young's modulus for the material. 
In seeking better refractory materials for furnace walls, we require to 
keep this product as low as possible. A low Young's modulus generally 
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means a correspondingly low tensile strength, but if a material could be 
found with a lower modulus, and a strength not reduced, it would, other 
things being equal, be superior. The search for a material with a low 
dilatation with rise of temperature seems more promising, and it is probably 
on account of its low dilatation that silica brick has proved so useful in 
furnace construction. 



The Two-dimensional Slow Motion of Viscous Fluids, 

By L. Bairstow, F.E.S., Professor of Aerodynamics, Imperial College of 
Science and Technology, Miss B. M. Cave, and Miss E. D. Lang, M.A. 

(Eeceived June 30, 1921.) 

The present paper is a contribution to the treatment of problems which 
require a solution of the differential equation v^-^/r = 0. Amongst such 
problems are to be found not only the very slow motions of a viscous fluid in 
two dimensions, but also the flexure of thin flat plates.^" 

The prosecution of the investigation has been made possible by the support 
of the Department of Scientific and Industrial Eesearch, which has provided 
financial assistance to enable two of us, to devote the whole of our time to 
the research, and our thanks are offered to the Department for its assistance. 
We also desire to acknowledge the facilities afforded by the Governing Body 
of the Imperial College of Science and Technology in placing a room at our 
disposal in the Department of Aeronautics. 

The method of attack was suggested by the results of an earlier paperf on 
the solution of Laplace's equation, y^i/r = 0. It was there found that, for 
any forms of single or double boundary, the solution of problems on the 
torsion of cylinders or the irrotational motion of an inviscid fluid could be 
made to depend on a number of definite integrals. In general, the evalua- 
tion of these integrals involves the use of graphical and mechanical methods, 
and much progress has been made in the direction of simplification of the 
processes. It is hoped that a paper may be presented elsewhere in the near 
future showing the application to a number of engineering problems. 

The solution of the equation v^'^ == given in the present paper follows 
the earlier work in its generality as to boundary forms and also in the 
general dependence on graphical and mechanical integration for arithmetical 

■^ Lamb, ' Hydrodynamics,' p. 604. 

+ ' Eoy. Soc. Proc.,' A, vol. 95, pp. 457-475 (1919). 



